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Abstract 

The spatial non-locallity (dispersion) of the transport equations results in a nonli- 
near dependence of the voltage drop U on the distance between the points of measuring. 
Therefore the results of the usual two-probe measurements of the conductivity depend 
essentially on the relation between the sample linear size L and the spatial disper- 
sion scale R of the generalized diffusion coefficient D(q,u). This makes it possible to 
get information on the character of the spatial non-locality of D(q,ui) in the Ander- 
son localization regime and, in particular, on the value of the correlation multifractal 
dimension D2 of the electron wave functions near the mobility edge. 

PACS 71.30.+h, 72.15.Rn 

1. In the last few years considerable attention has been focussed on the 
problem of spatial dispersion (non-locality) of the kinetic coefficients of dis- 
ordered systems in the Anderson localization regime [|TJ— ^j] . The cause of 
such interest was the realization of the fact that the character of the qu- 
dependence of the generalized diffusion coefficient D(q : cu) near the mobility 
edge is intimately connected to critical behavior of the electron wave functions 
and, in the end, is determined by the scenario of the metal-insulator transi- 
tion []I|,|3|. Indeed, the Berezinskij - Gorkov localization criterion |7j requires 
that in the localized phase D(q, 0) vanish simultaneously for all values of the 
wave number and, according to the one-parameter scaling hypothesis, the 
relation 

D(t; q, uj) = b 2 - d D{b 1/v t] bq, b d cu) , 2 < d < 4 (1) 

must be satisfied ||. Here t = (S — S c )/S c is the distance to the mobility 
edge S c , 6 is a scaling factor, d is the space dimension, v is the correlation 
length critical index. 

There are two qualitatively different versions of the critical behavior of 
D(q,cu) that obey these general requirements. 
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• According to Chalker's hypothesis [jlj], the multifractal structure of the 
electron wave functions near the Anderson transition (t — > and/or to — > 0) 
leads to an anomalously strong spatial dispersion of the generalized diffusion 
coefficient, whose scale is R = min(^,L CJ ) — > oc, where £ oc \t\~ v is the 



correlation length, L u oc ^D{uj)/uj oc u; -1 ^ is the electron diffusion length 
in a time of ~ 1/uj Depending on the relation between £, and g 

Chalker [jlj] distinguishes four main types of the asymptotic behavior: 

D( q ,u) = D (i) d ~ 2 { {qR) h-^ g*J; fl = ndn(f,L u ). (2) 

Here Do is the Drudian diffusion coefficient, / is the mean free path, 77 is a 
critical index related to the correlation multifractal dimension of the wave 
functions D 2 (r] = d- D 2 ) 0. 

• On the other hand, Suslov's symmetry approach to the localization the- 
ory H predicts supression of the spatial dispersion of the diffusion coefficient 
in the vicinity of the Anderson transition down to atomic scales ~ Xp. More 
recently, this conclusion was confirmed within the framework of a general- 
ized formulation |5],|6j] of the self-consistent Vollhardt - - Woffle theory [|IT| . 
According to 0,0], in the Anderson localization regime 

%«) = i 7p, iR«h 0) 



where the non-locality scale R oc y/Dor\D(t, u)/Dq\ and decreases as D(t,u) (x 
[min(£, L ul )] 2 ~ d until saturation is reached at R ~ A^. 

In the survey by Suslov |3j] it is pointed out that the absence of anomalously 
strong spatial dispersion of the generalized diffusion coefficient near the mo- 
bility edge does not contradict the concept of multifractality of the electron 
wave functions, it only indicate that the equality 77 = d— 2 (or D 2 = 2) should 
be satisfied. The well-known Wegner's result 77 = 2e (e = d — 2 <C 1) 0] is 
directly connected with the critical behavior of the inverse participation ra- 
tio. At the same time, the relation between this quantity and D(q,u) used 
in [P, |2], pD|1 cannot be considered correct for several reasons.^ Therefore, in 
our opinion, the above contradiction is only apparent. The same may be 
said about the results of numerical modelling 77 = 1.2 ± 0.15, 77 = 1.3 ± 0.2, 
77 = 1.5 ± 0.3 j2j and 77 = 1.3 ± 0.2 jjn| obtained by different methods for 



d = 3. 

This dilemma, touching upon the fundamental notions of the Ander- 
son localization, calls for both theoretical and experimental solution. In 

lr The detailed discussion of the basic arguments pro and con hypothesis 77 = d — 2 may be found in S. 
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this letter we derive a material equation relating the current density in the 
spatially nonuniform case to the experimentally measured difference of the 
electrochemical potentials, and propose a measuring scheme that allows one 
to obtain information on the degree of non-locality of the diffusion coefficient 
of charge carriers. 

2. The voltage drop U measured in the spatially nonuniform case is equal 
to the difference in electrochemical potential AU = A(/e = Acp + A/i/e 
between the corresponding points of a conductor. Therefore in calculating the 
current density it is necessary to take into account the response of the system 
to both the mechanical perturbation (electrical potential ip) and the thermal 
one (chemical potential fi) caused by the nonuniform electron distribution 
induced in the conductor.^ 

Applying the general equations of the linear response theory to 



the problem under consideration, we obtain a material equation relating the 
current density to the electrochemical potential gradient, for q <C kp its 
Fourier representation has the form 

j{Q,v) = -iqen F D(q,u)({q,cu) = -iqa{q,u)U{q,u) j 'e , (4) 

where np is the density of states at the Fermi level, a(q,uj) is the measured 
electrical conductivity connected with D(q : u) by the Einstein relation. It 
should be emphasized that a(q, to) is not Kubo's kinetic coefficient 

a(q, u) = e n F , 5 

1 + i{q z /uj)D{q,uj) 

which, unlike a(q,u) (see Eq. (g)), relates the total current density to the 
electrical field acting in the system. It is in the uniform case (q = 0) only 
that a(0,u) = a(0,u), and Eq. (g) coincides with the Einstein relation. 

Eqs. (Q) and (|5j) describe the non-local linear response of a spatially 
unbounded and homogeneous system. In the general case, the non-locality 
of material equations is of a more complicated character. In the typical 
experimental situation, the sample has the shape of a plane-parallel slab of 
thickness L. In this case the diffusion propagator of charge carriers G(x, x'\ to) 
is a solution of the equation 

d q 
— iuG(x, x'\ lS) — — J D(x, y\ u)—G(y, x'; co)dy = 5(x — x') (6) 

ox_ L/2 dy 



2 The nonuniform spatial electron distribution results not only in perturbation of the electrical field in 
the conductor, which is involved in if, but also in the appearance of the diffusion term in the measured total 
current density. Just the latter term is the response of the system to the thermal perturbation. 
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with the open boundary conditions 

G(x,x';u)\ x ^ =±L / 2 = 0. (7) 

The integral kernel of this equation (non-local diffusion coefficient) enters into 
the material equation that relates the electrical current to the electrochemical 
potential gradient (voltage drop) 

L/ r 2 ~ i 9 
I(x) = —e 2 npS I D(x : x';u)——U(x')dx' : (8) 

-L/2 QX 

where S oc L is the area of the sample cross section. The solution of the 
boundary problem (§D, (0) m the absence of spatial dispersion (D(x, x'\ uS) = 
D(u)5(x — x') ) is given by well-known method of images \\T4\ 



G(x,x';lu) = [G(x — x + 2nL\ uS) — G(x + x + (2n + 1)L; to)} , 

n=— oo 

n( \ 1 exp(iqx) 

G{x]uj) = y E ■ ~ 2n , n • 9) 
L q —ilu + q l D[u) 

Here G(x, ui) is the diffusion propagator of charge carriers in a spatially 
unbounded homogeneous system. We assume that the solution (0) holds true 
in the presence of spatial dispersion if its scale is small in comparison with the 
sample size (i?<Ci). In this case the integral kernel of the material equation 
(§) is expressed, by analogy with (0), in terms of the generalized diffusion 
coefficient D(q,u) of the spatially unbounded and homogeneous system 

~ 2 +°° 

D(x, x'\ to) = — E [D{q n , oj) sm.q n x sin q n x'+D(q n ,uj) cos q n x cosq n x'} , (10) 

where q n = 2im/L and q n = n(2n + 1)/L are discrete values of the wave 
numbers. 

Taking into account that the voltage drop U(x) is an odd function of x, 
and the current strength I(x) = I is constant along the studied sample, we 
will represent them through the corresponding Fourier series. Then, upon 
substitution of (|TD|) in the material equation (g), it is not difficult to find the 
expression for the voltage drop Fourier coefficient U n . Finally the x-depen- 
dence of the voltage drop can be presented as the following Fourier series: 

4/ +oo f_l\n+l 

U ( X ^ = T^~2 — ^ ^7v Tsmq n x, \x\ < L/2 . (11) 
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3. When the diffusion coefficient is independent of the wave number q n 
the series ([TTp gives, for a finite-size conductor, the usual conductance def- 
inition g(L,u) = L d_2 e 2 n j p J D(0, u) = L d ~ 2 a(u) and describes the linear in- 
dependence of the voltage drop U(x) within the sample (\x\ < L/2). The 
spatial non-locality of the generalized diffusion coefficient D(q n ,u) changes 
the conductance definition as: 

1 8 +°° 1 
= y (12) 

g(L,u) LSe 2 n F n = qlD(q n ,u;) 

and results in a nonlinear x-dependence of the voltage drop flTTl). Thus, 
information on the spatial dispersion of the generalized diffusion coemcient 
of charge carriers D(q : cu) can be obtained by measuring the nonlinear part 
of U(x) (O): 

AU(x) = U(x) - ^U(L/2) . (13) 

Lj 

Consider a sample in the form of a plane-parallel slab of thickness L ^> R 
(R is the spatial non-locality scale) with ideal ohmic contacts on the opposite 
surfaces and with two potential measuring probes which symmetrically lo- 
cated at a distance x from the mid-plane (see Fig. [5(a)). It is best to situate 
the measuring probes near the x max points where AU (x) reaches its maximum 
value AU max . Another pair of potential probes, the role of which is played 
here by the ideal ohmic contacts, are necessary for measuring the non-linear 
part of the voltage drop (p~3|) . 

Substituting (|3]) into ([IT]), ( [TBI ) yields AU(x) = 0. Strictly speaking, this 
equality is fulfilled in the case of ideal compensation of the linear part of 
the voltage drop. In any case, however, the measured signal is small in pa- 
rameter R/L <C 1. The nonzero contribution to AU(x) comes only from 
the nontrivial g-dependence of the generalized diffusion coemcient from (@). 
Fig. [I|(b) shows typical dependences of AU (x) calculated using linear inter- 
polation between Chalker's asymptotic formulas (|2|) for the inverse diffusion 
coemcient 1/D(q,u) at qR <C 1 and qR ^> 1. In this case x max does not 
depend on the scale of spatial non-locallity of D{q 1 uS) and assumes values on 
the interval x max « (0.70 -r- 0.76)L/2 for 77 = 1.1 1.5. The non-linear part 
A^max calculated at these points has the following asymptotic behaviors 

ac/ — (!)"-{ ;T r -, £>!; L -i« L - ^ 

Two ways for measuring the signal may be suggested. The first is 
to analyze the frequency dependence of At/ max (o;) in a sample with a fixed 
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Figure 1: (a) The scheme of the measurement of the nonlinear part of the voltage drop 
(here is depicted a longitudinal section of investigated sample); (b) Typical dependence of 
the nonlinear part of the voltage drop on the potential probe position x. Curve is calculated 
using Chalker's asymptotic formulas (0) for 77 = 1.3 and for R/L = 0.1. 



disorder level in a small enough vicinity of the mobility edge on the metallic 
side of transition (\t\ <C 1). According to and the predictions of Ref. 
as the frequency decreases, an increase in AC/ max oc uj~ 
be observed until saturation is reached at A£7 max oc 



■n/d 



[LU > UJr 



(f/i) 



V 



) should 

<C u c ). 



The typical frequency cu c is determined by the relation L UJc = £ or Ku c « 
X F d n F l \t\ vd , where Xp is de Broglie's wave length at the Fermi level. 

The second way is to measure the dependence of A£7 max ( |Ilf) on a dimension- 
less distance t to theimobility edge S c at a fixed frequency to. The stress- 
tuning technique [|T5| seems to be the simplest method of changing t in the 
vicinity of S c . A very suitable material for such measurements is Si:P, in 
which this technique makes it possible to attain \t\ ~ [|T5|. Accord- 

ing to (0) and the predictions of Ref. 0, as t decreases, an increase in 
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AC/max oc t~ vr] (t ^> tj) should be observed until saturation is reached at 
A[/ max oc (L^/L) 71 (t <C tu). Here the typical value t w is determined by the 
relation L w = £ or XpUphu « I^T^- 

Such dependences obtained by Eqs. QTTQ, (HH) using the interpolation 
R 1 = + Lj 1 for the non-locality scale, are plotted in Fig. 0(a,b). Since 
in a sufficiently small vicinity of the mobility edge the non-locality scale of 
D(q,cu) (|2|) takes on anomalously large values R = min(L w ,^), the value of 
AC/ max (HI) is quite accessible for measurement. The estimates show that for 
Si:P samples with the typical phosphorus concentrations np ~ 10 18 cm~ 3 and 
for really attainable values \t\ « 10~ 2 10 -3 , similar dependences of AU max 
(see Eq. ( |T^ ) and Fig. g) should be observed in the frequency region accessi- 
ble for probe measurements. For example, at t ~ 10~ 3 the typical frequence 
value is uj c « 10 3 10 4 s _1 (see Fig. @(a)), and the corresponding correlation 
lengths are of order ~ £ ~ 10~ 3 ^- 10" cm. 

So, we believe that the existence or supression j5|-|5f of the anomalous 
spatial dispersion of the generalized diffusion coefficient near the Anderson 
transition can be verified experimentally. 

The authors are grateful to I.M. Suslov for having drawn their attention 
to this problem and for helpful discussions. 
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